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METHOD FOR DETERMINING TJ IE^ ylAGE ^QUAL JTY 
OF AN OPTICAL IMAGING SYSTEM 

[0001] For example, a closed first diffraction ring at the edge of the first Rayleigh region 
can be considered as a sign of diffraction-limited optics. It is disadvantageous that this 
assessment must be considered only an integral information. And more specific quantitative 
information about the distribution of the rest of the imaging errors to different types of error, 
such as a sphere, coma or astigmatism, can also not be gained in this way. 

[0002] In another procedure, the fit of the individual optical components is checked by 
interferometry to arrive at assertions concerning the geometric errors, e.g., of a lens body, 
which can then be converted to system-dependent imaging errors. 

[0003] In this connection, system-dependent influencing factors are also already detected 
insofar as the measurement wavelength of the interferometer conforms to the working 
wavelength or wavelength spectrum of the illumination light. In more complicated optical 
systems, specially adapted interferometers are also occasionally used to check the image 
quality under given constraints and at the correct working wavelength with respect to the total 
imaging system. 

[0004] This is applied, for example, in imaging objectives for steppers or scanners to be 
used in semiconductor microlithography. This procedure requires a relatively high technical 
complexity and is therefore very cost-intensive and not usually employed in connection with 
microscope manufacture. 

[0005] Further, it is known to measure the wavefront of optical imaging systems with so- 
called Hartmann or Shack-Hartmann wavefront sensors or with sensors operating on similar 
principles. This also calls for relatively elaborate technology and, for that reason, the 
corresponding measuring systems are usually only designed to provide measurements only 
for different subsystems which, however, have similar interfaces, e.g., microscope objectives 
for microscopy. 

[0006] In this connection, a continual problem, particularly for microscope producers, is 
the lack of an available general testing procedure that would make it possible to determine 



the image quality as accurately as possible for the different optical imaging systems which 
differ from one another with respect to optical, geometric and mechanical parameters. 

[0007] Also, this problem exists not only during the manufacture and adjustment 
processes, but also in connection with quality control of imaging systems that are already in 
use by the customer. 

[0008] Further, the determination of the image quality for a plurality of field positions of 
the imaging system is uneconomical or inaccurate in all of the previously known procedures. 

[0009] The following source literature is cited in this connection: Joseph Geary, 
"Wavefront sensors", SPIE Press 1995, and Daniel Malacara, "Optical Shop Testing", Wiley 
Verlag 1992. 

[0010] Proceeding from this prior art, it is the object of the invention to further develop 
the previously known methods in such a way that a more accurate quantitative determination 
of the image quality of optical imaging systems is possible in a more economical manner and 
the results of the determination of image quality can also be used to determine the influence 
of initially unknown samples on the amplitude distribution and phase front distribution of the 
illumination light, of which the amplitude distribution, in particular, is known. 

[001 1] According to the invention, the following method steps are provided in a method 
of the type mentioned above: 

adjusting the subassemblies relative to one another in such a way that it is 
possible to image a sample on the detection device; 

recording a plurality of images of the sample from different adjusting or 
reference planes near the focus plane, wherein the detection device is adjusted relative to the 
image plane, the sample is adjusted relative to the object plane, or the objective is adjusted 
relative to the sample; 

improving the image quality by means of image processing, particularly to 
reduce noise, to compensate for local variations in sensitivity of the detection device, and to 
center the intensity centroids respectively on a predetermined location in the images; 

computational linking of spatially resolved image information, of adjustment 
values and system variables relating to the optical imaging system, and of information 
concerning the sample with the aim of determining characteristic numbers that are 
characteristic of the wavefront deformation caused by the imaging system; and 



outputting the characteristic numbers and associating them with the imaging 
system for description of the image quality. 

[0012] The basic idea of the invention consists in initially making a plurality of 
recordings, a so-called image stack, of a sample, preferably a sample with a pinhole, with the 
operative optical imaging system, of which imaging parameters such as wavelength, 
numerical aperture, sigma (as degree of coherence of the illumination), scanning, image size 
(in pixels), magnification, and so forth, are known, wherein the individual recordings of this 
image stack come from different focusing distances, i.e., from different defocusing 
adjustments near the focus plane. 

[0013] In other words, the sample is imaged repeatedly, each time with different 
defocusing, on the spatially resolving detection device and the image information that is 
obtained by pixel is stored. This provides a stack of layer images one on top of the other in 
which every image is made up of a multitude of image information in the form of intensity 
values. 

[0014] The individual images of the image stack are initially subjected to measures for 
improving the image quality which are known, per se, in the technical field of image 
processing. These measures relate particularly to 

reducing noise or improving the signal-to-noise ratio; 
compensating for local sensitivity differences of the individual sensor 
elements (pixels) of the detection device; 

linearizing the dynamic range of the sensor elements; 
subtracting background noise; 

centering the intensity centroids in the individual images on a predetermined 
location; and/or 

cutting out relevant image portions, e.g., from axial points and field points, 
while preventing edge cropping. 

[0015] These measures for improving image quality are particularly advantageous with 
respect to the accuracy that can be achieved with the additional method steps according to the 
invention. 

[0016] After improving the image quality, the image information and all relevant data 
such as adjusting values and system variables relating to the optical imaging system and 
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information concerning the sample (type and characteristics of sample) are transferred to the 
evaluating device for combining based on predetermined computing programs. 

[0017] Along with the transfer to the evaluating device, it is decided at the same time 
whether or not to carry out a deconvolution of the image information in order to separate the 
sample influences from the influences of the imaging system and illumination. This decision 
is made depending, for example, on the diameter of a pinhole d PH in the sample compared to 
the object-side Airy diameter d A iry of the imaging system and on the accuracy to be achieved. 
The corresponding control parameters for the deconvolution, such as pinhole diameter or, 
e.g., parameters for a Tikhonov regularization, are likewise entered into the evaluating device 
in case it is decided to carry out the deconvolution. 

[0018] Further, the following must be determined prior to the evaluation or provided for 
the evaluation: 

the desired accuracy of the evaluation or of the results of the evaluation; 

the quantity of characteristic numbers which are to be determined in the 
assessment of the image quality or upon which the assessment is to be based; 

the selection of a computing program to be used for the evaluation from a 
stock of available computing programs; 

termination criteria for the evaluation, for example, the magnitude of the 
residual error, the number of iteration cycles, or the duration of the evaluation; 

whether the evaluation is to be carried out in one step or in multiple steps; in 
the latter case, the aim is to arrive at results faster, to increase convergence and stabilize the 
process; and 

whether or not the defocusing in the object space, i.e., the adjustment of the 
sample relative to the object plane or the adjustment of the objective relative to the sample, 
must be taken into account during the evaluation; in this case, e.g., the parameter space for an 
iteration calculation to be provided must be increased so that influences originating from the 
object space defocusing can be taken into account. 

[0019] In a further development of the operating sequence, pupil functions can be taken 
into account, e.g., the results of the measurement of the real illumination by pupil imaging 
using the Bertrand system or the theoretically determinable apodization of the pupil as a 
characteristic of the design of the imaging system. 
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[0020] In a particularly preferred embodiment of the invention, the characteristic numbers 
are determined by analytic evaluation in a first step and by iterative evaluation in a 
subsequent second step, and the results from the first step are transferred to the second step as 
starting values and are further processed until a given termination criterion is reached. 

[002 1 ] The first step is used for quickly finding the target area of the results for the 
characteristic numbers with a robust method, but can already deliver results with limited 
accuracy. The second step, which is based on the, as yet, inexact results of the first step or 
which makes use of these results as starting values, serves to determine more exact results. 

[0022] In the first, analytic step, a first set of characteristic numbers which are equivalent 
for the image quality are determined, e.g., in the form of Zernike polynomials up to a desired 
order of magnitude. For this purpose, the information or data introduced in the evaluation are 
analyzed, for example, by the extended Zernike method. 

[0023] Generally, the characteristic numbers resulting from the analytic evaluation do not 
yet satisfy the requirements for accuracy in the assessment of the image quality of an optical 
imaging system. Therefore, the analytic determination is followed by an iterative process. 

[0024] The iterative evaluating process substantially incorporates the image information, 
the image recording parameters of the optical imaging system, the characteristic numbers 
determined in the preceding analytic evaluation step for image quality, and residuals as 
residual errors between retrieved and measured intensity distribution. 

[0025] The iterative phase of the evaluation can be carried out either pixel by pixel with 
error minimization, e.g., based on merit functions or derived method variants, or two- 
dimensionally (taking into account the image contents). The iterative method to be used is 
preferably based on the Gerchberg or Gerchberg-Saxton principle, optionally associated with 
the Levenberg-Marquardt, Gaussian, or non-least-square optimization methods, known in 
technical circles. 

[0026] In addition, it is also conceivable that the iterative evaluation is not preceded by an 
analytic evaluation but, rather, that the image information and specifications of the imaging 
system are furnished directly to the iterative evaluation, wherein the image information which 
is not dependent upon the sample and is possibly deconvolved is taken into account. 
However, starting with the analytic evaluation has the considerable advantage that starting 



values are available for iteration which appreciably increase the convergence and the 
certainty of actually finding the absolute minimum of residuals. 

[0027] When defocusing is carried out during the recording of the image stack in the 
object space and the object space defocusing is also to be taken into account in the 
evaluation, the parameter space for the iteration must be increased, preferably doubled in 
linear development, so that one or more parameters describing the defocusing can be 
associated with every characteristic number for the image quality and the defocusing effects 
are accordingly taken into account in the iteration. In the linear case, the characteristic 
numbers then have the following form: 

Ci = c if + Az*c i)d , 

where c i>f is the characteristic number of the image quality, Az is the advance in direction of 
the optical axis, and Ci,d is the defocusing term. 

[0028] The data and information upon which the iteration is to be based can be taken over 
selectively from different antecedent steps of the determination of measurement values or 
from the results of theoretical system calculation. 

[0029] For a robust and reliable or stable evaluation, known global optimization methods 
can be incorporated in the iterative evaluation, e.g., simulated annealing or self-learning 
computing processes. 

[0030] The output of the evaluation results is provided in the form of characteristic 
numbers that describe the image quality, preferably as Zernike coefficients, as image 
information, in the form of residuals as residual errors between retrieved and measured 
intensity distribution, and as defocusing parameters associated with the determined Zernike 
coefficients. 

[003 1] Further, the method according to the invention is used for calculating the imaging 
errors that are now also known along with the image quality from the images of any, initially 
still unknown, samples recorded by the same imaging system. 

[0032] For this purpose, the image information obtained with the initially still unknown 
samples is subjected to post-processing in which the characteristics of the imaging system are 
taken into account by deconvolution. In this way, the specific device characteristics are taken 
into account and also corrected at the same time when imaging the sample. It is 
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advantageous when deconvolving the characteristics of the imaging system from the sample 
images, to correct the influence of the finite object size from the point images at the same 
time. 

[0033] In this way, it is also possible to factor in the influence, e.g., of a stepper in 
microlithography in the characteristics of an image by convolution again in order to correct 
the characteristics of the measuring system and accordingly obtain an image such as would 
have been generated by the stepper. 

[0034] The change of reference plane for the recording of the image stack should 
advantageously be carried out in steps to be determined. 

[0035] The image stack can be recorded by means of an individual spatially resolving 
detection device, wherein an adjustment of the detection device or of the sample is required. 
Alternatively, it is also possible to record the image stack using a plurality of spatially 
resolving detection devices which must then be arranged in the desired reference planes. In 
the latter case, it is possible to record all images of the image stack simultaneously and also, 
depending on the control, to record the images belonging to an image stack successively in 
time. 

[0036] When the image quality is to be determined in relation to different positions in the 
visual field of the imaging system, it can be provided in a development of the invention that 
the determination of the characteristic numbers is based on a plurality of samples that are 
positioned adjacent to one another in the sample holder. Accordingly, for each recording, the 
samples supply image information that can be associated with the respective positions. 

[0037] It is likewise possible to arrange a sample with a plurality of objects in the visual 
field of the imaging system so that information concerning the image errors in the visual field 
can be determined in an image stack. In arranging the sample with respect to the 
characteristics of the imaging system, it is advantageous to ensure that the objects do not 
influence one another in the arrangement. 

[0038] The samples can be binary objects, i.e., pure amplitude objects, e.g., lithography 
masks. 



[0039] In another particularly preferred embodiment of the invention, the determination of 
all of the data and information supplied in the evaluation and the evaluation, including the 
output of the evaluation results, take place automatically. 

[0040] In some cases, it is advantageous when obtaining the image stack to use an 
exposure control which, in spite of the changed reference planes, ensures an optimal 
illumination of the sample so that the signal-to-noise ratio in the images is optimized already 
during the recording. 

[0041] For illumination of the sample, it is also conceivable to use a laser beam having a 
beam waist in the object plane, which typically leads to a low sigma value and also to a 
Gaussian intensity distribution in the pupil of the imaging system. 

[0042] In a further development of the method according to the invention, the 
determination of the influence of initially unknown samples on the amplitude distribution and 
phase front distribution of the illumination light, of which this distribution is known, is 
provided with the following method steps: 

determining the wavefront deformation characterizing the optical imaging 
system based on a sample whose characteristics are known; 

exchanging this sample for a sample which is to be examined and whose 
characteristics are initially still unknown; 

determining the wavefront deformation again while taking into account the 
influence of the sample to be examined after the method steps described above; 

separating the characteristic numbers for the image quality resulting under the 
influence of the unknown sample and the characteristic numbers for the image quality that 
were determined without the influence of the unknown sample; and 

determining the characteristics of the initially unknown sample from the 
results of this separation. 

[0043] Determining the characteristics of samples in the manner mentioned above is 
connected to the preceding method steps and assumes an operative imaging system for which 
the characteristic numbers characterizing the image quality are already available. With 
respect to this imaging system, the image parameters, such as wavelength of the illumination 
light, numerical aperture, sigma, scanning, magnification, image size (in pixels) and 



defocusing parameters (associated with the characteristic numbers for the image quality), are 
also known or can be determined. 

[0044] An, as yet, unknown sample to be examined is placed in this imaging system with 
known imaging characteristics, or an image of the sample is projected therein. A series of 
images, that is, an image stack, of the sample to be examined is recorded. In this case also, 
the focus distance is changed for each image as was already described. A defocusing- 
dependent exposure control can also be used for this purpose to adapt the illumination of the 
sample to the respective focus plane. 

[0045] The images of the sample to be examined which are obtained in this way are 
processed with respect to their quality and the signal-to-noise ratio is again improved, e.g., by 
dark image correction, the pixels are linearized (flat-fielding), the background noise is 
subtracted, the intensity centroid is centered in every image, relevant image portions are cut 
out and/or noise filtering can be provided, so that an image stack with information-containing 
images that can be evaluated is supplied for further processing. 

[0046] The centering of the intensity centroid using the data values contained in the image 
can be carried out directly, or by means of correlation since the images do not now terminate 
at the edge of the camera field (in contrast to the recording of an image stack for determining 
the characteristic numbers of the pure imaging system), or can be carried out purely 
mechanically when the lateral sequence in the X,Y-direction has been determined beforehand 
in relation to the z-position for the imaging system in the respective reference plane, possibly 
also with interpolation of intermediate planes. 

[0047] The image information resulting after this processing, the information concerning 
the imaging system, the information concerning the adjustment with respect to each image, 
and the characteristic numbers for the image quality, are now transferred to an evaluating 
device. 

[0048] The following are supplied, in particular, for the evaluation: 

the desired accuracy of the evaluation or of the evaluation results; 

the quantity of characteristic number for describing the image quality, e.g., the 
quantity of Zemike coefficients upon which the evaluation and results are to be based; 

selection of the computing methods and optimization methods to be used 
insofar as the evaluating device has a plurality of suitable computational algorithms; 
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preset of the termination criterion for an iterative evaluation, e.g., residual 
errors, number of iterative cycles, duration of evaluation; 

the one-step or multiple-step evaluation; the aim of the multiple-step 
evaluation is to shorten the computing time, increase the convergence and the certainty of 
finding results in the iteration, or also to stabilize the computing process; 

taking into account an object space defocusing that may have been carried out, 
for which purpose the parameter space for the iteration calculation should be increased in 
order to be able to take into account a variation of the characteristic numbers of the image 
quality that result from the object defocusing; and 

a criterion for the transition from an analytic evaluation to an iterative 
evaluation, e.g., residual errors, duration of evaluation, and so on. 

[0049] According to the invention, the characteristic numbers for the image quality of the 
imaging system are again determined analytically with the, as yet, unknown sample in a first 
step of the multiple-step evaluation. For this purpose, Zernike polynomials are determined 
up to a desired order, e.g., according to the extended Zernike method. 

[0050] The iterative determination of the characteristic numbers for the image quality is 
carried out in a subsequent second step of the evaluation, wherein the results of the preceding 
analytic evaluation are used as starting values for the iteration in order to increase the 
convergence of the iteration method in this way. 

[0051] Input quantities for the iterative step of the evaluation are again the image 
information, the characteristic numbers for image quality (e.g., in the form of Zernike 
coefficients) and, for purposes of comparison, residuals as residual errors between retrieved 
and measured intensity distribution in the images of the sample, and, as the case may be, the 
system parameters for the defocusing characteristics. 

[0052] The iteration can again be carried out either pixel by pixel with error minimization, 
e.g., according to the Gerchberg-Saxton method, followed by optimization methods, e.g., the 
non-least-square method, Levenberg-Marquardt, or the like. 

[0053] For a robust and reliable or stable evaluation, known global optimization methods 
can be incorporated in the iterative evaluation, e.g., simulated annealing or self-learning 
computing processes. 



[0054] It is preferable to use iterative methods operating pixel by pixel because 
information must be obtained about a sample that is initially still unknown. The use of the 
two-dimensional iterating method is recommended preferably with samples having only 
extensive structures and relatively small gradients, e.g., when examining certain biological 
samples. 

[0055] In an embodiment of the method according to the invention, the object space 
defocusing is taken into account and, in connection with this, the parameter space for the 
iteration is enlarged so that every characteristic number determined for the image quality can 
be associated with defocusing parameters and the defocusing effects are taken into account in 
the iteration. The values of the imaging system can be used by themselves as starting values. 

[0056] The evaluation results in characteristic numbers of the image quality of the 
imaging system which are influenced by the sample to be examined and which consequently 
deviate from the characteristic numbers that were determined beforehand with a known 
sample. 

[0057] The following procedure can be carried out for separating the sample 
characteristics from the characteristics of the imaging system: 

[0058] 1. A best focus determination is carried out for the obtained image stack and, 
based on the characteristic numbers of the image quality of the imaging system, a simulated 
focus series of very simple object structures, e.g., of a point image, an edge, or column, is 
calculated only for the imaging system at the focus positions of the measured image stack. 
By deconvolving the images of the measured image stack with corresponding images of the 
simulated image stack, the images of the measured image stack are freed from the 
characteristics of the imaging system . The characteristic numbers of the image quality for 
the sample by itself are determined from this measured image stack. 

[0059] With the sample-oriented characteristic numbers of the image quality determined 
in this way, it is possible to calculate backward to the amplitude distribution, phase 
distribution and/or intensity distribution in the object plane for the sample by itself in order to 
obtain the sample characteristics. In principle, a measured image series can also be used 
instead of the simulated image stack, although computation will be less accurate. 

[0060] 2. The characteristic numbers of the imaging system, as were determined with 
the known sample according to the method steps described in the beginning, are subtracted 
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from the characteristic numbers of the imaging system that were determined under the 
influence of the sample to be examined. In this way, the sample characteristics are separated 
and can be outputted in suitable form, e.g., in the form of information concerning spatial 
amplitude distribution or intensity distribution, spatial phase distribution, as isoamplitudes, 
isointensity surfaces, or isophase surfaces. 

[0061] Particularly when using pixel-oriented calculation methods, an unwrapping of the 
phase surfaces may be necessary. In this case also, it is possible to calculate backward to the 
amplitude distribution, phase distribution or intensity distribution in the object plane with the 
sample-oriented characteristic numbers determined in this way in order to make the sample 
characteristics more precise. 

[0062] 3. Based on the characteristic numbers of the image quality for the imaging 
system with the known sample on the one hand and the characteristic numbers of the image 
quality of the imaging system without the sample to be examined on the other hand, the 
optical field distributions in the object plane are determined and deconvolved by calculating 
backward and the information concerning the characteristics of the sample by itself are 
obtained. 

[0063] The output of results is preferably carried out in the form of image information, 
characteristic numbers for image quality, e.g., Zemike coefficients, defocusing parameters for 
the measured Zernike coefficients, information concerning amplitude, phase or intensity that 
can be calculated back to the influence of the sample. 

[0064] In principle, the output of results is carried out in the same way in the three 
variants indicated above, but different additional information can be outputted at the same 
time depending on the selected variant. 

[0065] The method according to the invention will be described in the following with 
reference to an embodiment example. In the accompanying drawings: 

[0066] Fig. 1 shows a schematic view of optical subassemblies used in an optical imaging 
system; 

[0067] Fig. 2 shows a simplified view of an image stack comprising a plurality of images 
of a sample obtained from different focus settings; 
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[0068] Fig. 3 shows an example for the recording of an image stack with defocusing 
within the image space; 

[0069] Fig. 4 shows an example for the recording of an image stack with defocusing 
within the object space; 

[0070] Fig. 5 shows a view of a practical optical imaging system; 

[0071] Fig. 6 illustrates the basic procedure for evaluation by iteration; and 

[0072] Fig. 7 shows proof of the image quality of an optical imaging system according to 
Fig. 5 in the form of a value table. 

[0073] The subassemblies according to Fig. 1 substantially comprise an illumination 
device 1, a sample holder 2 with a sample 2.1 which is positioned in or near the object plane 
3, an objective 4, and a CCD camera 5, as spatially resolving detection device, which is 
positioned in or near the image plane 6 of the objective 4. The subassemblies have a 
common optical axis 7 which does not necessarily extend in a straight line. 

[0074] The subassemblies are adjusted relative to one another in a first method step in 
such a way that it is possible for the sample 2.1, which can exist physically or in the form of 
an image, to be projected onto the receiving surface of the CCD camera 5. The receiving 
surface of the CCD camera 5 comprises an array of sensor elements (pixels), and information 
concerning the intensity of the impinging illumination light influenced by the sample can be 
read off at the output thereof. 

[0075] In a second method step, images are recorded from different defocusing ranges and 
result, in their entirety, in an image stack as is shown schematically in Fig. 2. The individual 
images are designated by A, B, C, D and E. The image stack extends in the direction of the 
optical axis 7 over a depth corresponding to the sum of distances d a b to d de . 

[0076] Of course, the quantity of images is not restricted to the quantity shown here; 
rather, substantially any quantity can be chosen. However, it is recommended that an odd- 
number quantity be selected in order to obtain, by specific adjustment, an image from the best 
focus plane and the rest of the images from defocusing ranges with defocusing symmetric to 
the best focus plane. 

[0077] Each of the images A to E comprises a grid of image information arranged in rows 
i and columns j. This image information is made up of intensity values, whose magnitude 
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corresponds to the output signal of an associated sensor element (pixels) of the receiving 
surface of the CCD camera 5, and represents the image of the sample. 

[0078] In Fig. 2, the image information Ay, By, Cy, Dy and Ey which is arranged 
successively in direction of the optical axis 7 ideally represents one and the same area of the 
sample recorded from different focus planes and, therefore, also with different intensity 
values. 

[0079] The direction of the optical axis 7 corresponds to the coordinate direction Z of the, 
ideally, Cartesian coordinate system, while each of the images lies in the plane defined by the 
coordinates X and Y. 

[0080] In principle, the image information can be obtained by means of a VIS microscope, 
a UV microscope or another imaging system. 

[0081] As is shown in Fig. 3, the recording of the image stack can be carried out, for 
example, by adjusting the CCD camera 5 relative to the image plane 6, i.e., the recordings are 
obtained with defocusing within the image space. The distance between sample 2.1 and 
objective 4 remains unchanged. Depending on the different reference planes, a pinhole in the 
sample 2.1 causes the intensity caustic shown in Fig. 3. 

[0082] Alternatively, the defocusing can also be achieved by changing the distance 
between the sample 2.1 and the objective 4, that is, by defocusing in the object space as is 
shown in Fig. 4. In this case, the distance between the receiving surface of the CCD camera 
5 and the objective 4 or pupil plane 8, for example, remains constant. 

[0083] The adjustment of the sample in direction of the optical axis 7 and the selection of 
different defocusing positions result in different wavefront deformations in the pupil plane 8 
which lead to different intensity sections that are associated with the individual defocusing 
positions as is shown in Fig. 4. 

[0084] In the next method step, the quality of the images and the image information is 
improved, for example, according to known methods of image processing, in that the noise is 
suppressed, only the intensity values above a threshold are taken into account for the 
subsequent evaluation, and measurement data are filtered and/or subjected to equalization. 

[0085] Accordingly, local sensitivity differences of the CCD camera 5 are compensated at 
the same time. Further, the intensity centroid of every image is centered on a predetermined 
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point in the image so that the intensity centra ids lie on a straight line, for example, in 
positions Ay, By, Cy, Dy and Ey (see Fig. 2). 

[0086] In the following, the invention will be explained in more detail with reference to an 
optical imaging system, according to Fig. 5, comprising the subassemblies which are shown 
schematically in Fig. 1. 

[0087] The following is assumed for this optical imaging system: 

the sample has a pinhole with a diameter d PH = 300 nm; 
illumination light with the wavelength of 248 nm is used; 
the pixel size at the sample is 45 nm; 
the numerical aperture of the imaging system is 0.2; 

the illumination aperture corresponds to the numerical aperture of the imaging 

system; 

the illumination of the sample is carried out with partially coherent light at 

sigma- 0.8; 

the diameter of the Airy disk in the image is 1 .5 12 jam; 
the depth of focus is 6.2 jxm; 

the defocusing from image to image is carried out within the depth of focus 
range at ±1 RE (RE = Rayleigh unit), ±3 RE and ±0.8 RE or ±6.2 jim, ±18.6 pm and ±5 pm; 
and 

the selected number of images is 2 1 . 

[0088] These data and the image information from the image stack are taken as the basis 
for the subsequent evaluation in which the following procedure is followed in practice: 

the characteristic numbers for the image quality of the imaging system are 
determined analytically by the extended Zernike method; 

two-dimensional iteration is carried out using Marquardt-Levenberg 
optimization in iteration steps, shown, for example, in Fig. 6; 

object space defocusing is partially taken into account. 

[0089] The results of the evaluation are indicated in the value table in Fig. 7 which also 
shows, among other things, the static Zernike values depending on an indexing after fringe 
scaling 1 to 25 as a description of the image quality. 
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[0090] This presents concrete evidence of the image quality for the imaging system 
according to Fig. 5. This evidence can be supplied to the customer, e.g., as a certificate, upon 
delivery of a corresponding device for further use or can be made available for customer 
service. For example, the quality of the imaging system can be determined and documented 
also on location at regular intervals or on demand. 

[0091] The method according to the invention has been explained primarily in connection 
with microscope imaging systems. Of course, it is also possible and advantageous to apply 
the method according to the invention in connection with other systems, particularly with all 
systems which image at a finite distance, such as binoculars, projectors, geodetic optical 
devices, camera systems, photographic apparatus, medical observation instruments and, as 
the case may be, also in stepper optics and scanner optics. In principle, it is possible to apply 
the method according to the invention for all optical systems that image a sample or an 
object. The method according to the invention is extensively independent of wavelength and 
can therefore be used in principle for all wavelengths. 

[0092] The expression "wavefront deformation" used in the description of the invention 
relates to deviations relative to a predetermined wavefront; it is not always desirable that an 
imaging system generates a plane wavefront. By extension, this also applies to the 
expression "imaging error". 

[0093] The possibility of graduated method steps according to the invention is particularly 
advantageous for determining the characteristic numbers in analytic and iterative steps. 
Further, graduation of the iterative method steps, per se, is also conceivable in that the 
iteration is carried out in multiple steps. 

[0094] As regards allowing for object space defocusing in the determination of results, it 
should be noted that the iteration suggested according to the invention should be carried out 
initially without additional parameters for the object space defocusing so as to arrive at the 
range of desired result values as quickly as possible. When the range of resulting values is 
found, e.g., characterized by transition parameters such as the ratio of residuals to signal-to- 
noise ratio in the respective image, quantity of iteration loops, duration of iteration, slight 
deviations in the results of an actual iteration cycle relative to the preceding iteration cycle, 
etc., the computing method should be expanded by parameters relating to the object space 
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defocusing, and the parameter space of the variables can also be enlarged at the same time, or 
doubled in case of linear development of the object space defocusing. 

[0095] Directly presetting the object space defocusing can also disadvantageously lead to 
errors because of the many parameters to be preset. For this reason also, iteration can 
initially be carried out without taking into account the object space defocusing, and the object 
space defocusing is only taken into account in later, likewise iterative, method steps. 
Improved stability and convergence of the computing process can be achieved in this way. 

[0096] Further, it is also conceivable to directly allow for the object space defocusing 
analytically, although this can also lead to the possible error rate, or to allow for the object 
space defocusing analytically followed by one or more iterative steps and, depending on the 
preset, initially without allowing for the object space defocusing thereafter. 

[0097] The image quality need not necessarily be represented in the form of Zernike 
coefficients as is indicated by way of example in the description of the invention. When the 
image quality is described in the form of Zernike coefficients, there is no information 
contained therein concerning, e.g., scattered light, contrast, distortions, image sequence in 
defocusing, transmission, field, image field curvature, image inclination, or dispersion. 
Insofar as this information is also significant or at least desirable with respect to the 
assessment of imaging systems, it can be determined by further examination based on the 
measured and simulated image stacks of known samples or from the image centering. 



Reference Numbers 



1 illumination device 

2 sample holder 

3 object plane 

4 objective 

5 CCD camera 

6 image plane 

7 optical axis 

8 pupil plane 



